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ABSTRACT 

This paper is the third in a series on tests of gravity using observations of stars and 
nearby dwarf galaxies. We carry out four distinct tests using pubhshed data on the 
kinematics and morphology of dwarf galaxies, motivated by the theoretical work of 
Hui et al. (2009) and Jain & Vanderplas (2011). In a wide class of gravity theories 
a scalar field couples to matter and provides an attractive fifth force. Due to their 
different self-gravity, stars and gas may respond differently to the scalar force leading 
to several observable deviations from standard gravity. HI gas, red giant stars and main 
sequence stars can be displaced relative to each other, and the stellar disk can display 
warps or asymmetric rotation curves aligned with external potential gradients. To 
distinguish the effects of modified gravity from standard astrophysical phenomena, we 
use a control sample of galaxies that are expected to be screened from the fifth force. In 
all cases we find no significant deviation from the null hypothesis of general relativity. 
The limits obtained from dwarf galaxies are not yet competitive with the limits from 
cepheids obtained in our first paper, but can be improved to probe regions of parameter 
space that are inaccessible using other tests. We discuss how our methodology can be 
applied to new radio and optical observations of nearby galaxies. 

Key words: Modified gravity 



1 INTRODUCTION 

The explanation of the observed accelerated expansion of 
the universe remains a mystery, but possible solutions to this 
puzzle fall into two broad categories. The first is to posit the 
existence of a new component of the universe with an appro- 
priate equation of state to cause the observed acceleration. 
This "dark energy" may be due to an exotic particle or field, 
or be related to the vacuum energy of space itself. A second 
approach is to seek to explain this acceleration through mod- 
ifications to the field equations of general relativity (GR) 
itself. In recent years this modified gravity (MG) approach 
has received attention and various approaches have been ac- 
tively developed - for a review see e.g. Jain & Khoury (2010). 
Even on theoretical grounds, modified theories of gravity are 
well motivated: GR is unlikely to be the complete theory of 
gravity owing to its singularities and non-renormalizeability. 
Hence it is generally considered to be the low energy effec- 
tive action of some UV-complete theory, though our interest 
is modifications in the long distance/low energy regime. 

A modification of GR on large (astrophysical) scales 
generically leads to scalar-tensor theories of gravity, where 
a new scalar field couples to gravity. Equivalently these the- 



ories can be described via a coupling of the scalar field to 
matter, which leads to enhancements of the gravitational 
force. Nonrelativistic matter - such as the stars, gas, and 
dust in galaxies - will feel this enhanced force, which in 
general leads to larger dynamically inferred masses. The dis- 
crepancy can be up to a factor of 1 /3 in f (R) or DGP gravity. 
We note that photons, being relativistic, do not feel the en- 
hanced force, so that lensing measurements probe the true 
mass distribution. 

This enhanced gravitational force should be detectable 
through fifth force experiments, tests of the equivalence prin- 
ciple (if the scalar coupling to matter varied with the prop- 
erties of matter) or through the orbits of planets around 
the Sun (Will 2006). Khoury k Weltman (2004) proposed 
that nonlinear screening of the scalar field, called chameleon 
screening, can suppress the fifth force in high density envi- 
ronments such as the Milky Way, so that Solar System and 
lab tests can then be satisfied. This screening was originally 
suggested to hide the effects of a quintessence-like scalar 
that forms the dark energy and may couple to matter - 
generically such a coupling is expected unless forbidden by 
a symmetry - so there are theoretical reasons to expect a 
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screening effect to operate in either a dark energy or modi- 
fied gravity scenario. Indeed, since there are only a handful 
of screening mechanisms (Vainshtein and symmetron screen- 
ing are discussed below) , small scale tests of gravity that rely 
on distinct signatures of screening are useful discriminators 
of cosmological models. 

In this paper we will consider deviations from GR exhib- 
ited in theories that rely on chameleon screening. Qualita- 
tively similar behavior occurs in symmetron screening (Hin- 
terbichler & Khoury 2010) and the environmentally depen- 
dent dilaton (Brax et al. 2010) and the tests we present here 
apply to these mechanisms as well^ . 

While local dynamical studies suggest that any fifth- 
force effects must be screened within the Milky Way, 
dwarf galaxies in low-density environments may remain un- 
screened as the Newtonian potential (which determines the 
level of screening) can be an order of magnitude smaller 
in magnitude than in the Milky Way. Hence dwarf galaxies 
can exhibit manifestations of modified forces in their infall 
motions, internal dynamics and segregation of stars, black 
holes and gas. Hui et al. (2009), Hui & Nicolis (2010), Jain & 
Vanderplas (2011), Chang k Hui (2011), Davis et al. (2012) 
and Jain et al. (2012) discuss various observational conse- 
quences. Indeed the small scale dynamics of nearby galaxies 
may have a bigger signal of MG than large-scale perturba- 
tions (see Jain 2011, for a discussion of observational ap- 
proaches) . 

In this work, we discuss observable effects on disk galax- 
ies that arise due to their internal dynamics or their inter- 
action with a neighbor or other sources of an external grav- 
itational field. We focus on late-type dwarf galaxies as these 
are most likely to be unscreened. The observational effects 
are rooted in the fact that compact objects like stars within 
the dwarf galaxies can self-screen, and thus respond to a 
different gravitational potential than the diffuse unscreened 
components like gas and dark matter. So, for example, if the 
stellar component of a dwarf galaxy self-screens, it will lag 
the dark matter and HI gas disk in the infall towards another 
galaxy. This may lead to a separation of the stellar disk from 
the center of mass of the dark matter and from the HI disk, 
and lead to observable distortions of the morphology and/or 
dynamics of the stellar disk. Jain & Vanderplas (2011) com- 
puted a variety of dynamical consequences of these effects. 

However, the origin of modified gravity effects may be 
difficult to disentangle from other astrophysical processes. 
To address this, we create a control samples of screened 
galaxies which are not expected to show any of the ex- 
pected modified gravity effects. By comparing statistics of 
the unscreened sample and screened control sample, we may 
derive constraints on modified gravity parameters. The di- 
vision of observed galaxies into screened and unscreened 
samples is accomplished based on an estimate of the local 
value of the external gravitational potential at the location 
of each galaxy: see Cabre et al. (2012) for a description of 
the methodology used to distinguish the screened and un- 



^ In order to identify the screened and unscreened galaxies it is 
necessary to calibrate on N-body simulations, which give differ- 
ent criteria for different screening mechanisms. While we do not 
implement criteria for any mechanism other than chameleon, our 
results appear robust to the details of screening. 



screened samples. Jain et al. (2012) used this methodology 
and obtained modified gravity constraints based on cepheid 
variable stars and other distance indicators. This paper is 
the third in the series of papers on astrophysical tests of 
gravity: 

In this paper we perform several of the proposed tests of 
modified gravity using currently available data. Because of 
the low signal-to-noise of these types of tests, a large amount 
of high quality data will be required to obtain competitive 
constraints on modified gravity models. Therefore, one of 
the main goals of this paper is to show the strengths and 
weaknesses of currently available data for these tests, and 
explore the requirements for future data in this area. 

This paper is organized as follows. In section 2 we de- 
scribe the method to generate potential map of the local 
Universe, which allows us to prepare our samples of screened 
and unscreened galaxies. In section 3 we briefiy summarize 
each of the observable effects we will use to test modified 
gravity. In section 4 we describe the offset between stellar 
and gaseous components, which includes: offset between HI 
and optical centroids (section 4.1), offset between optical 
center and Ha center (section 4.2), and offset between red 
giants and main sequence stars (section 4.3). In section 5 we 
study MG-induced warping of edge-on galaxies. In section 
6 we describe the tests based on gaseous velocity rotation 
curves, which includes HI - Ha velocity differences (section 
6.1), and asymmetry in the Ha curves (section 6.2). We con- 
clude in section 7. 



2 POTENTIAL MAP OF THE LOCAL 
UNIVERSE 

Each of the modified gravity tests presented in this paper re- 
lies on comparing various observationally-derived statistics 
of a screened and unscreened sample of galaxies. Therefore, 
the first step of the analysis is to determine whether each 
galaxy in our sample is likely to be screened, i.e. the 'fifth' 
force due to the scalar field is suppressed by the deep po- 
tential well of the galaxy or its environment. Because the 
screened sample should be free of any fifth-force effects, we 
use it as a control sample. The detailed procedure for deter- 
mining screening has been described and tested via N-body 
simulations of f{R) gravity in Cabre et al. (2012); we sum- 
marize these results here. 

For an isolated spherical halo in f{R) gravity, the 
chameleon effect is suppressed if fno ^ 2/3 |(/)Ar|/c^, where 
(/)Ar is the Newtonian potential of the object (Hu & Sawicki 
2007) . We will approximate galaxies as isolated spherical ha- 
los, and make use of this condition to determine whether an 
object is self-screened or environmentally screened: i.e. we 
look for the galaxies which satisfy the following conditions: 

l^intl 3 |0ext| 3 

> ^JRO or > -ino (1) 

The first condition describes the self- screening condition, 
with |(/)int| = CM/rvir, where the M and rvir are the mass 
and virial radius. We typically use measurements of the cir- 
cular velocity at large radii to establish self-screening, since 
we work with rotating disk galaxies for our tests. The sec- 
ond condition describes environmental screening condition, 
with |(/)ext| representing the Newtonian potential due to the 
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neighbor galaxies within the background Compton length 
Ac- The external potential |(/)ext| is evaluated using neigh- 
boring galaxies as: 

di<Xc+ri 

where di is the physical distance to the neighboring galaxy 
with mass Mi and virial radius . Ac is related to the back- 
ground field value fno (Schmidt et al. 2009): 

Ac « 32ySMpc (3) 

Identifying screened and unscreened regions requires 
knowledge of the masses and three-dimensional positions of 
all the objects in the sky. Since our knowledge is limited 
by finite survey sizes and depths, we assume that most of 
the contribution to the environmental screening potential 
comes from clusters and groups of galaxies. This is a rea- 
sonable assumption as the potential of a rich cluster is up 
to two orders of magnitudes larger than a typical galaxy. 
Cabre et al. (2012) compiled a catalog of bright clusters 
with known redshifts based on Abell (Abell et al. 1989) and 
ROSAT (Ebeling et al. 1996) observations. We remove many 
clusters from the Abell catalog as they lack spectroscopic 
redshift or have redshift greater than 0.1. Our final cata- 
log contains 675 clusters of galaxies with galactic latitude 
|6| > 20. 

The next most significant contribution comes from 
galaxy groups. Yang et al. (2007) provides a nearly com- 
plete catalog of groups in the SDSS region, including around 
400,000 galaxies covering the redshift range 0.01 — 0.2. In 
addition, we use the group catalog from the 2M++ all sky 
survey (Lavaux & Hudson 2011) to create a combined cat- 
alog of 153380 SDSS and 3984 2M++ galaxy groups. The 
final contribution to the catalog are local galaxies within 10 
Mpc, identified by Karachentsev et al. (2004). 

With this catalog in place, an estimate of the mass of 
each object is required in order to calculate the potential. We 
determine the mass using either the mass-luminosity relation 
given by Reiprich & Bohringer (2002) or the mass- velocity 
dispersion relation given by Evrard et al. (2008). The masses 
of SDSS groups and local galaxies are taken directly from 
Yang et al. (2007) and Karachentsev et al. (2004), respec- 
tively^ . 

Using the mass of the clusters, groups and galaxies, we 
generate a potential map of the sky. In this paper we use 
values of the Compton length between Ac = 3.2 Mpc and 
Ac = 1-0 Mpc corresponding to fno parameters 10~^ and 
10~^. Due to survey incompleteness, the derived potential 
map is only an approximation for many regions of the sky. 
This could potentially lead to systematic inaccuracies in la- 
beling of screened and unscreened objects. A detailed dis- 
cussion of these uncertainties can be found in Cabre et al. 
(2012). 



^ Note that Karachentsev et al. (2004) reports the 

mass for only 313 of the 451 galaxy sample. See 

http:/ /www. sas.upenn.edu/~vinu/screening/index. html for 
an extensive description of the catalog. 



3 OBSERVABLE EFFECTS MODIFIED 
GRAVITY 

In chameleon gravity theories, objects with large gravita- 
tional potential may be self-screened (see Eq. 1). This may 
lead to observable differences in the dynamics of different 
galactic components: the gaseous disk, dark matter halo, and 
giant stars may feel the fifth force, while the self-screened 
main sequence stars move according to standard gravita- 
tional dynamics. The different forces on these components 
lead to several potentially observable morphological and 
kinematical differences between the two components. In this 
section, we summarize the expected effects of the fifth force, 
as well as potential observational tracers of these effects. 

3.1 Summary of physical effects of modified 
gravity 

The following effects of modified gravity are expected to ap- 
pear in unscreened dwarf galaxies. Quantitative estimates 
of the magnitude of these observables are given in the fol- 
lowing sections. For more details and motivation of these 
predictions, see Jain & Vanderplas (2011). 

(i) Off*set between stellar and gaseous components 
(section 4): Chameleon theories predict that in the pres- 
ence of a suitable potential gradient, low density, unscreened 
components like dark matter or neutral hydrogen clouds 
(HI) will fall faster along an external potential than will self- 
screened stars. This difference in acceleration may cause a 
measurable offset between the centroid of the stellar disk 
and the gaseous disk. The population of red giant stars, be- 
ing partially screened, may also show an observable offset 
from the population of main sequence stars. The size of the 
offset depends strongly on the mass and concentration of 
the dwarf galaxy's halo, as well as the distance to neighbor- 
ing galaxies. Jain & Vanderplas (2011) estimate this offset 
may be several percent of the virial radius of suitable dwarf 
galaxies. We carry out three tests of such offsets in section 
4. 

(ii) Warping of the stellar disk (section 5): As the 

halo moves along an external potential gradient, it pulls at 
the lagging stellar component. If the potential gradient is 
aligned with the axis of rotation, the steep slope of the halo 
potential may introduce a U-shaped warp in the stellar disk. 
The shape and strength of the warp depends on the mass 
profile of the dwarf galaxy halo, as well as the strength of the 
external potential gradient. The warp is expected to align 
with this potential gradient. Jain & Vanderplas (2011) es- 
timate the warp to be of order 0.1 kpc, though it varies 
significantly depending on the external fifth force acting on 
the dwarf galaxy. 

(iii) Rotation curve tests (section 6): If a disk galaxy 
falls edge-on, or more generally if the external potential gra- 
dient is not aligned with the galaxy's axis of rotation, the 
offset between stellar and halo centers can perturb the stellar 
disk and cause an asymmetry in the stellar rotation curve. 
This is due to the fact that the dominant force on the stellar 
disk is from the potential of the halo, rather than the poten- 
tial of the disk itself. The asymmetry can be as high as 10 
km/s. The asymmetry can be distinguished from standard 
astrophysical sources by comparing the stellar disk with the 
gaseous disk. A version of this test is show in section 6.2. 
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In principle the rotation velocity of an unscreened gaseous 
disk is also enhanced relative to the stellar disk by a factor of 
Y^l + AG/G, where AC is the contribution of the fifth force 
to the self-gravity of a galaxy (Hui et al. 2009). We find that 
current data is unsuitable to carry out this test: see section 
6.1 for a discussion of current data and its limitations. 

In isolated cases, any of the above observable properties 
could potentially be explained via standard astrophysical ef- 
fects. But under modified gravity we would expect to see 
multiple effects in conjunction, and to see offsets and warps 
aligned with the direction of the external gravitational force. 
More importantly, we expect these effects to be more preva- 
lent in unscreened galaxies than in screened galaxies: such a 
correlation would be difficult to explain using conventional 
explanations, such as bars, spiral arms, asymmetric drift, 
and external tidal effects. In the following sections, there- 
fore, we explore aggregate statistics relating to the presence 
of these observable effects in screened and unscreened pop- 
ulations. 



3.2 Observational Tracers 

In order to trace gaseous components of dwarf galaxies, we 
use observations of the 21 cm line which results from the hy- 
perfine transition of neutral hydrogen. In order to trace the 
stellar component, we use the distribution of optical light. 
Traditionally, Ho; is also be used as a kinematical tracer of 
the disk, but this presents some complications for studies of 
modified gravity. 

Ho; emission in dwarf galaxies arises because the high 
energy radiation from stars ionizes the neutral gas in the 
ISM, forming a "Stromgren Sphere" . As the ions within this 
region recombine, the photons emitted in the n = 3 n = 
2 transition of Hydrogen lead to the Ha emission line. 

Historically, this line has been used to measure the stel- 
lar rotation curve, and is expected to give accurate mea- 
surements of the dynamics of the stellar disk in the absence 
of modified gravity. However, the possibility of a fifth force 
which affects stars and gas calls this expectation into ques- 
tion. Recall that the star is screened by its own potential, 
but this potential may not be strong enough to screen the 
whole of the surrounding sphere. 

A rough calculation shows that for a star of mass 
IOM0 and radius 5Rq, the surface potential is ~ 4 x 10~^. 
Therefore, this star will be screened for background field 
value less than 4 x 10~^. For fm = 1 x 10~^ (2 x 10~^), 
this potential is only sufficient to screen ionized gas within 
4.6 X 10"^ (2.3 X 10"^) parsec of the star - that is, about 4 
(20) times the star's radius. For typical HII regions which 
have radii of order lOpc, the region screened by the stellar 
potential is negligible. Therefore, we would expect that the 
HII region feels the fifth force and orbits faster than the as- 
sociated star - this raises interesting questions about what 
fraction of stars remain associated with significant HII re- 
gions. Thus, under the assumptions of modified gravity, it is 
not safe to assume that Ha traces the dynamics of the stel- 
lar component. We will simply assume that Ha is a gaseous 
tracer and leave a detailed study for later work, though we 
provide a brief empirical test Section 6.1. 

A better tracer of the rotation of the stellar component 
under modified gravity may be molecular absorption lines 



in the stellar spectra, which arise much closer to the star. 
Absorption lines such as Mgib or Call would trace the stel- 
lar rotation curve accurately under both modified gravity 
and general relativity. However, these lines are difficult to 
observe, particularly in the low-surface-brightness outskirts 
of the galaxy. 



4 OFFSET BETWEEN STELLAR AND 
GASEOUS COMPONENTS 

In this section we consider three tests based on the first ob- 
servable effect: the expected offset between unscreened and 
self-screened components of the unscreened galaxy sample. 
Simple calculations indicate that under reasonable approxi- 
mations, this offset should be large enough to be observable 
(Jain & Vanderplas 2011). We consider three potential av- 
enues toward this observation below. 



4.1 Offset between HI and optical centroids 

We first compare the optical centroid, which traces the stel- 
lar component, to the centroid of HI gas measured via its 
21cm emission. Here the optical observations are taken from 
the SDSS r-band (Abazajian et al. 2009) while the HI cen- 
troids are taken from the ALFALFA survey which used the 
Arecibo radio telescope (Giovanelli et al. 2005, 2007; Sain- 
tonge et al. 2008; Kent et al. 2008; Martin et al. 2009). The 
astrometric precision (about 24 arcseconds) of these observa- 
tions are limited by the relatively large beam size of the radio 
maps. Though the low resolution of the HI image makes any 
offset difficult to measure for a single image, the statistical 
uncertainty can be reduced by considering a large sample of 
galaxies. 

To reduce the uncertainty in the spatial centroiding, 
we limit our sample to those galaxies within 40 Mpc. To 
minimize the error introduced by the incompleteness in our 
catalog described in Section 2, we further restrict our analy- 
sis to the SDSS region, where we have information on galaxy 
groups down to a mass of 5 x lO^^M©. Finally, to limit the 
effects of different detection thresholds between the SDSS 
and Arecibo images, we keep only the galaxies whose opti- 
cal and HI centroids are within 1 arcmin. After these cuts 
are made, the resulting catalog contains 245 dwarf galaxies. 

In order to identify galaxies which are large enough to 
be self-screened, we use the maximum velocity estimated 
from W50, i.e. the velocity width of the source line pro- 
file measured at the 50% level of each of the two peaks. 
This was corrected for instrumental broadening by the AL- 
FALFA team. However, this must be further corrected for 
turbulence and inclination. The turbulence correction was 
performed according to Tully & Fouque (1985) using an av- 
erage velocity dispersion in the gas component of 8 km/s 
(Begum et al. 2006; Geha et al. 2006). Finally, we apply the 
inclination correction based on the observed axis ratio q in 
the SDSS r-band. The maximum velocity is found via: 

W50c ... 

Vuia.^ = , , (4) 

2^1^ 

where VF50c is the turbulence corrected W50 and q is the 
axis ratio of the galaxy. For a background field value fno, 
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Figure 1. The distribution of the measured offset between the HI and optical centroids for galaxies observed by the ALFALFA and 
SDSS surveys. The black (red) points represent screened (unscreened) galaxies for model parameter values /rq = 1 x 10~^ (left panel) 
and /ro = 2 X 10~^ (right panel) 



the unscreened galaxies satisfy (see Jain et al. 2012) 

V 100 7 ^ 2 X 10-7" ^ ^ 

The above criteria yields 68 (85) unscreened and 177 (160) 
screened galaxies for /rq = 2 x 10 ~^ {/rq = 10 ~^). 

In Figure 1 we show the measured Hi/optical offset for 
samples of different screening level, along with the distri- 
butions of offsets for two values of /rq. The errors shown 
are estimated from the 24 arcsec centroid error of the radio 
maps, converted to a physical distance using the distance to 
each galaxy. 

The interpretation of any observed offset is made diffi- 
cult by the fact that there are several potential astrophysi- 
cal (non-MG) sources of separation, as well as the fact that 
we cannot predict the magnitude or direction of the MG- 
induced separation. We address this by proposing an aggre- 
gate statistical model of the data, which takes into account 
the expected difference between the screened and unscreened 
samples. 

To account for non-MG related separation, we assume 
that standard astrophysical effects will produce a scatter in 
the distribution of projected Hl-optical separations centered 
at zero and with a standard deviation acn- This scatter will 
be present in both the screened and unscreened samples. To 
account for MG-related separation, we assume that within 
the unscreened sample, there will be an additional Gaus- 
sian scatter centered at zero with standard deviation ctmg- 
Then the observed separation for a screened galaxy and an 
unscreened galaxy is respectively 



^unscr 
Si 



= abs{si^ 
abs{si^ 



(6) 



where sf^^ and s'^^^^^ are the observed offsets for the 



are (unobserved) values drawn from Gaussian distributions: 



^ M{Q,(T%R + (jIig) 

e. -A/'(0,a,') 



(7) 



where J\f denotes a Gaussian distribution and ai is the 
(known) measurement error associated with observation i. 

With this model in place, we can find the most likely pa- 
rameter combination within a Bayesian formalism, by writ- 
ing 



p[cFMG,CFGR\{Si}) p[cFmg)p[cFGr) 



n 



g{si]^,aGR + al) 



n 



^(s^;0,CrMG + ^GR + 



(8) 



Si 



'M{^,(jIr) 



where g{x;/j,a^) is the Gaussian probability distribution 
with mean /i and width a; p{(tmg) and p{(tgr) are the pri- 
ors on the parameters; and S^^^ and S^^^^^ are respectively 
the sets of screened and unscreened observations. 

In Figure 2 we show the 1 — , 2—, and 3 — a contours of 
the posterior distribution given in Eq. 8 as a function of ctmg 
and acR- The left panel shows the result for the observed 
data: this data is not precise enough to rule-out a ~ Ikpc 
scatter due to modified gravity: this is primarily due to the 
large observational error in the measurement of the radio 
centroid. In the right panel, we see that if the aggregate 
astrometric accuracy were to be improved by a factor of 
10 (i.e. a 2.4 arcsecond centroid accuracy, or a factor of 100 
increase in the sample size, or some combination of these two 
improvements), the observations would be sensitive enough 
to detect a scatter due to modified gravity of ^ 0.5kpc, 
on the order of the offsets predicted by (Jain & Vanderplas 
2011) for typical MG theories. 

One important detail of these results for such low signal- 
to-noise data is their dependence on the precise determi- 
nation of the measurement uncertainty. Within the model 
given by Equation 8, the measurement error is degenerate 
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Figure 2. The 1-, 2-, and 3-<t contours for the scatter in observed Hi/optical centroid separations in nearby dwarf galaxies. The left 
panel shows the results given the observed sample of 245 dwarf galaxies described in §4.1. Because of the large beam size of the radio 
maps, the data cannot distinguish between modified gravity and general relativity. The right panel shows a simulated sample drawn 
from the model with a scatter due to modified gravity of 0.5 kpc (indicated by the 'X'), with a factor of ten improvement in astrometric 
precision (or a combination of improved astrometry and larger sample size). The resulting posterior can distinguish the predicted effects 
of modified gravity at the l-cr level. 



with the astrophysical offset: any increase in assumed mea- 
surement error will be offset by a commensurate decrease in 
the astrophysical scatter at maximum likelihood. Further- 
more, because the uncertainty of the measurement domi- 
nates the signal by a factor of up to 4 for each galaxy, even 
a 10% overestimate or underestimate in the centroid error 
will lead to significant differences in the maximum likelihood 
fit. This caveat should be kept in mind when interpreting 
Figure 2. 

Thus, while the current data are not precise enough to 
either confirm or rule-out modified gravity theories, a large, 
combined radio / optical survey of nearby dwarf galaxies with 
very precise and well-characterized astrometry could provide 
strong constraints on typical modified gravity theories. 



4.2 Offset between optical center and Ho; 
kinematical center 

A primary weakness of the measurements explored in the 
previous section is the astrometric precision of centroid mea- 
surements in the radio. In this section we measure the same 
effect, but using the kinematical center of the Ha disk as a 
tracer of its centroid. Here we assume that the gas produc- 
ing the Hq! line is unscreened and traces the dark matter, 
which means that the rotation curve should (in principle) 
have the same center as the dark matter halo (see Sec 3.2). 
As above, we search for systematic offsets between the (self- 
screened) stellar component and this (unscreened) gaseous 
component. 

For this section, we use data presented by Persic & 
Salucci (1995): they re-analyze 967 Ha velocity rotation 
curves originally observed by Mathewson et al. (1992), and 
publish the offset between the kinematical center (which 
symmetrize the rotation curve) and the optical center. We 
divide these galaxies into a screened and unscreened sample 



using the methodology presented in Section 2, leading to 25 
(690) unscreened galaxies and 858 (193) screened galaxies 
for fRo = 2x 10-^ (fno = lO"^). 

In the left panel of Figure 3, we show the distance to 
each galaxy along with its mass indicator, the maximum ro- 
tation velocity. The mild correlation between these measure- 
ments can be traced to to observational bias: smaller, fainter 
galaxies are less likely to be observed at larger distances. 
Note that velocities range from 50 to 300 km/s, which allow 
us to test theories from /rq = 10~^ (galaxies self-screened 
if ^max > 225 km/s) to fno = 10"^ (^max > 70 km/s). The 
color of each point indicates whether the galaxy is screened 
at a level of fm = 10"^ ^ There is a relatively sharp bound- 
ary between the screened and unscreened samples seen at ^ 
120 km/s rotation velocity: this is due to the self-screening 
criterion met by larger-mass galaxies. 

Because the size of the expected MG displacement 
scales with the size of the dark matter halo (Jain & Van- 
derplas 2011), we work here with the relative offset S — 
Sohs/Topt where Sohs is the offset between the optical and 
kinematical centers and Vopt is the optical radius of the 
galaxy. This relative offset does not show an obvious cor- 
relation with galaxy rotation velocity, as can be seen in the 
right panel of Figure 3. 

To search for indications of modified gravity, we use 
the aggregate likelihood formalism developed in the previous 
section, replacing the absolute deviations s with the relative 
deviations 5. The resulting posterior is shown in Figure 4. In 
comparison to the ALFALFA results in Figure 1, the data 
posterior is much tighter: this is mainly due to the increased 
astrometric accuracy of these measurements. 

The left panel of Figure 4 shows results for the entire 
sample, and (for = 10~^) indicates an approximately 
1 — a detection of an offset due to modified gravity. Because 
this result uses the whole sample, however, it is potentially 



© 0000 RAS, MNRAS 000, 000-000 



Astrophysical Tests of Modified Gravity 7 



60 

20 


350 
300 
250 
200 
150 
100 
50 
0. 




-I 1 1 1— 




100 
75 
50 
25 


350 
300 
250 



Ng =671 
Nu =212 



20 40 60 80 100 
Distance (Mpc) 



100 
50 




120 140 0.00 



0.05 0.10 
Offset / Optical 



0.15 
radius 



0.20 



Figure 3. The distribution of galaxies used for the optical/kinematic offset test described in Sec 4.2. The left panel shows the distance 
and maximal rotation velocity (a tracer of mass) of the galaxy sample; the right panel shows the relative offset between the gaseous 
and stellar components. Black (red) points represent screened (unscreened) galaxies at /rq = 10~^'^. The clear mass dependence of the 
screening condition is due to self-screening of massive galaxies. 



biased by the sizes of the galaxies in the sample. As shown in 
Figure 3, the vast majority of screened galaxies in our sam- 
ple are not environmentally screened, but self-screened. The 
differential offset shown in this analysis, then, not only re- 
flects the difference between screened and unscreened galax- 
ies, but also the difference between large and small galaxies. 
Even absent modified gravity effects, we may expect statisti- 
cally significant differences in observed offsets between these 
samples. 

To limit the effect of this size bias, we repeat the anal- 
ysis, limiting our galaxy sample to those with a maximum 
rotation velocity less than 200 km/s (i.e. those galaxies be- 
low the dashed line in Figure 3. Having controlled for this 
bias, we see in the left panel of Figure 4 that the offsets are 
consistent with General Relativity for both screening levels. 

Still, this analysis does not rule out offsets due to mod- 
ified gravity with a scatter of approximately 2% of the op- 
tical radius of the galaxy. Under ideal conditions, Jain & 
Vanderplas (2011) showed that the expected offsets for typ- 
ical modified gravity theories could be as high as 10% the 
optical radius of the galaxy. This predicted offset decreases 
as the size of the galaxy increases, and as the distance to the 
neighboring galaxy increases. Thus a 2% upper-limit is not 
enough to definitively rule-out modified gravity theories. 



4.3 Offset between centroids of red giants and 
main sequence stars 

In modified gravity scenarios, red giant stars present an 
interesting case: their cores may be dense enough to self- 
screen, while their diffuse outer shells remain unscreened 
and feel the fifth force. This partial screening can lead to 
interesting effects: for example, Chang & Hui (2011) have 
recently shown that red giant stars evolve differently under 
modified gravity compared to main sequence stars, and that 



this difference can be used to constrain modified gravity the- 
ories. 

In this section we make use of red giant populations 
to constrain modified gravity in a different way: just as un- 
screened gas and dark matter will separate from the screened 
main sequence stars in the presence of an external potential 
gradient, partially screened red giants will be expected to 
separate from the population of main sequence stars as well. 

4.3.1 Data and centroid estimation of RGB and main 
sequence stars 

Recently, Dalcanton et al. (2009) published positions and 
magnitudes of individual stars within 69 nearby dwarf galax- 
ies (distance < 4 Mpc , |6| > 20 deg)^. Absolute magnitudes 
of these galaxies (Mb) range from -20 to -8. They observed 
nearly 10^ to 10^ stars in each galaxy. In our study we need 
full images of the galaxies to estimate the displacement be- 
tween RGB (red giant branch) and main sequence stellar 
populations as partial images of galaxies introduce artifi- 
cial displacement between those two populations. We use 
only those galaxies which pass this criterion. In addition, 
we exclude all the galaxies in which foreground contami- 
nation leads to poorly-defined red giant branches. Finally, 
we include galaxies only if they have at least 3000 observed 
stars. The final sample consists of 28 galaxies: under either 
screening assumption {/ro = 10~^ and /rq = 10~^) every 
galaxy is unscreened, which means this test cannot rely on 
a screened control sample. 

For each galaxy, we select red giants from the color- 
magnitude diagram. In Figure Bl we show an example of 
red giant star selection. Because the fifth force acts only on 
the outer regions of the largest giant stars, we tune the cutoff 
to select only the brightest, and therefore the largest, of the 

^ See http://www.nearbygalaxies.org 
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Figure 4. The Bayesian estimate of the contributions of MG and astrophysical effects to the offset between kinematical and optical 
centroids (see discussion of this model in Sec 4.1). The left panel shows a mildly significant offset due to modified gravity: we argue that 
this is due to the bias of the screened sample toward high masses. After controlling for this by removing galaxies with rotation greater 
than 200 km/s, the estimate is consistent with GR. 



red giant stars. Next we define the area of the galaxy within 
which the centroid will be calculated. This is done by taking 
isocontours of stellar number density significantly above the 
noise level. Errors in the centroids of the red giant and main 
sequence populations are calculated using bootstrap resam- 
pling with 50 realizations. These are shot-noise dominated. 
Since the number of red giants is at least an order of mag- 
nitude lower than the number of main sequence stars, the 
error in the centroid displacement is dominated by the error 
in the centroid position of the red giant sample. 



4-3.2 Estimation of the external fifth force 

Because of the lack of control sample, detecting modified 
gravity here depends on correlating the direction of displace- 
ment with the strength and direction of the background 
gravitational potential gradient at each dwarf galaxy due 
to its neighbors. In order to estimate this, we use the cata- 
log of nearby galaxies given by Karachentsev et al. (2004). 
This catalog is nearly complete out to a distance of 10 Mpc 
and it contains a total of 431 galaxies with known mass for 
313 galaxies. We fit a linear relation between the mass and 
magnitude (in logarithmic scale) to the galaxies with known 
mass and use this empirical relation to estimate the masses 
of the remaining galaxies. By using distances and masses of 
neighboring galaxies we can construct the force vector (per 
unit mass) as 



potential W can be calculated as 



^ GMj fj-r 
/-^ \ri—r\'^\ri—r\ 

zGunscr ' ' ' ' 



(9) 



where m is the studied galaxy mass at position r. Mi the 
mass of the neighbor z, which is at position fl and G is New- 
ton's constant. The sum includes only unscreened neighbors 
at a distance below the Compton length. The gravitational 



(10) 



where the sum includes all neighbors at a distance below 
the Compton length. Here we use the potential to estimate 
the magnitude and direction of the background gravita- 
tional force, as well as to estimate the potential to define 
the screening level for the dwarf galaxy sample. The force 
has some errors due to the uncertainty in distance (5%) and 
mass {(Tiog{M) — 0.2). In our case, uncertainty in mass is the 
dominant source of error. 



4.3.3 Results and comparison with theory 

In Figure 5 we plot the displacement of red giants from 
main sequence stars projected to the external force vector 
as a function of the force per unit mass. The force vector is 
previously projected to the plane of the sky perpendicular 
to the line of sight. In order to calculate the error in the pro- 
jected displacement, we take into account both the error in 
the centroids and the error in the vector force, and show the 
results for /^o = 10~^ and fno = 10~^. There is no statisti- 
cally significant correlation between the fore and projected 
displacement. The weighted average displacement in the di- 
rection of the external force for fno = 10~^ (fno = 10~^) 
is —10.2 ± 5.5 pc (—4.7 =b 6.4 pc), consistent with zero to 
within 2-cr (1-cr). The errors on these estimates include in- 
trinsic dispersion, as explained in Appendix A. 

In Figure 6 we compare the projected displacement to 
the prediction, given the external force and assuming a cored 
isothermal sphere (see Equation 2.3 of Jain & Vanderplas 
(2011) for the analytical expression of the offset). The pre- 
diction takes into account errors in the masses and distances 
of neighboring galaxies, as well as errors in the determina- 
tion of the galaxy mass, core radius, and core mass. The last 
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Figure 5. The displacement of red giants from main sequence stars projected along the external fifth force from a sample of 28 nearby 
dwarf galaxies. The displacement is plotted as a function of the external force per unit mass (Eq. 9). In the left (right) panel, we set 
Iro = 10~^(10~^). All the galaxies are unscreened. There is no significant correlation between the displacement and the applied force. 
The mean displacement in the direction of the external force is —10.2 zb 5.5 pc and — 4.7zb 6.4 pc for the two cases, consistent with GR. 



two quantities are estimated using the strong empirical cor- 
relation between virial mass and core radius and core mass 
(Swaters et al. 2003, 2011). We have assumed AG/G = 1 
for the predictions, but this result can be easily scaled for 
other values of AG/G. The mean difference between the 
observed displacement d and predicted displacement p is 
p-d= 12.5d=6.1 pc (-0.7±4.8 pc) for fno = 10"^ (10"^). 

We find that both these tests are unable to rule-out ei- 
ther GR or modified gravity. The main reason for this lack 
of discriminatory power is because most of the galaxies are 
very isolated, so the external force has a small magnitude. 
In Appendix B we discuss in detaila series of important sys- 
tematic uncertainties involved in this test. 



4-3.4 Future directions 

We find that tests of modified gravity through comparison 
of the positions of red giant and main sequence stars are 
inconclusive, mainly owing to the error associated with the 
small sample size. Also, because red giants are only partially 
screened, the signal may be diluted when compared to the 
tests of the previous sections. Additionally, the lack of a 
control sample means that the uncertainty in the external 
force is the limiting factor if one wants to constrain modified 
gravity with only the unscreened sample as such tests require 
much better prediction to compare the data. 

Future studies making use of this effect should focus 
on unscreened galaxies with very close neighbors, which will 
thus experience stronger external potential gradients, and 
drive the separation to more detectable levels. It would also 
be helpful to focus on similarly-sized screened galaxies to 
provide a control sample (though this pursuit may be limited 
by the relative paucity of nearby screened dwarf galaxies). 
The best approach would be to focus on a limited number 
of galaxies with regular morphology, to reduce intrinsic sys- 
tematic effects, with a very well defined external force. 

Another potentially interesting method for this sort of 
measurement would be to use integral field spectra of galax- 



ies and estimate the red giant and main sequence contri- 
butions to the spectrum as a function of position: any sys- 
tematic offset, aligned with the external fifth force in the 
inferred mass of red giant and main sequence stars could 
indicate the presence of a fifth force acting on the red giant 
population. 



5 WARPING OF EDGE-ON GALAXIES 

In the previous section, we discussed the expected offset 
of the stellar and dark matter components of suitable un- 
screened dwarf galaxies within modified gravity. As shown 
in Jain & Vanderplas (2011), this offset may lead to an ob- 
servable morphological signature: a bowl-shaped warp in the 
stellar disk, aligned with the direction of the external grav- 
itational force and this effect will be maximized when the 
gravitational force is perpendicular to the plane of the disk. 
When viewed edge-on, this will result in a symmetric U- 
shaped warp in the disk of the infalling galaxy. Because 
standard astrophysical dynamics (e.g. tides) could lead to 
similar effects, we use screened galaxies as a control sample 
and look for a statistically significant higher fraction of such 
U-shaped warps in the unscreened sample. 

5.1 Data and estimation of warp parameter 

We search for warps in low mass, nearly edge-on galaxies 
using the ALFALFA and SDSS data. The rotation velocity 
measured by the former is used to find self-screened galaxies, 
and images from the latter are used to estimate the galaxy 
warp. Because warp is difficult to observe in face-on galaxies, 
we limit our sample to those galaxies with an axis ratio of 
less than 0.6. In the resulting sample of 495 galaxies, 128 
(59) galaxies are screened and 367 (158) unscreened, for a 
Compton scale length of 3 (1) Mpc which corresponds to 
/^o = 1 X 10-^ (fno = 2 X 10-^). For /i^o = 2 x 10"^ we use 
only galaxies with Vmax < 200 km/s in the screened sample. 
In order to reduce the type of mass-related systematic 
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Figure 6. The estimated projected displacement, as in Fig. 5, is compared with the predicted displacement assuming AG/G = 1. In the 
left (right) panel, we set Jrq = 10~^(10~^). The solid line shows the one to one correspondence. It is evident that while the observations 
are consistent with GR, one cannot rule out modified gravity given the large errors in the predictions and systematics. 



effects discussed in the previous section, we include only 
smaller galaxies in which the self-screening condition (Eq. 5) 
does not apply: the separation of screened and unscreened 
galaxies is based only on the environment. The warp of the 
galaxies is found by slightly modifying the method used in 
Jimenez- Vicente et al. (1997): after sky subtraction, we first 
rotate the SDSS images so that the major axis of the galaxies 
align with the horizontal axis. The center of rotation is found 
using the SExtractor package (Bertin & Arnouts 1996). We 
mask stars and other contaminating nearby objects in the 
image. After removing all the contaminating objects and 
pixels in the frame we find the centroids for each column in 
the image by fitting a Gaussian to the column and taking the 
first moment. This leads to a 'warp curve', i.e. the centroid 
of the columns as a function of distance from the center 
of the galaxies. To clean the curve, we starting from the 
center of the galaxy and sweep through the curve, omitting 
any column which has a linearly-interpolated centroid that 
deviates more than three pixels from that of the previous 
column. An example of such a cleaned warp curve is shown 
on Figure 7. 

After estimating the warp curve, we quantify the overall 
warp by defining the following two measures: 



Wi — 



r 

JO 



x_ dx 
L L 



W2 — 



xydx 



- 2l r 



xydx 



(11) 



(12) 



for the left {x < 0) and right {x > 0) side of the galaxy. 

wi is dimensionless and proportional to the mean y/L 
displacement from the x-axis weighted by the distance from 
the center (x/L), where L = max(|x|). The normalizing fac- 
tor (l/L) assures that wi is independent of the choice of 
the measurement unit. W2 has units of kpc, and measures 
the average deviation from a horizontal plane of the galaxy 
weighted by x (horizontal axis). In each warp measure, the 



error is calculated via standard propagation of the width of 
each column centroid. 

We correct both warp measures for the inclination of the 
galaxy based on Geha et al. (2006). This is done by dividing 

the measured warp curve by sin(i) where i = 



for a 



relative disk thickness of 0.19. The warp is estimated within 
three half-light radii of the galaxy (beyond this radius, the 
images tend to become very noisy), though the results of 
this section do not vary appreciably as this limit is changed 
from 2 to 4 half-light radii. In Appendix G4 we define a 
third warp measure, and show that the following results are 
unchanged. 



5.2 Results 

We classify the warp as "U-shaped" if wi or W2 has oppo- 
site signs on the left and right sides, and "S-shaped" if the 
signs match. To quantify the warp strength, we average the 
absolute values of the warp in the left and right sides of 
the galaxy. Figure 8 shows the normalized distribution of 
wi (left panel) and W2 (right panel) for screened and un- 
screened galaxies displaying a U-shaped warp. The errors in 
the histograms are estimated using Poisson statistics in each 
bin. 

It is expected within modified gravity that unscreened 
galaxies with close neighbors will preferentially show large 
U-shaped warps. A K-S test shows that the wi distributions 
are similar between the screened and unscreened samples: 
this implies that the appearance of warp does not depend 
on the screening level. However, we find that W2 distribu- 
tions are marginally different: the unscreened galaxies are 
slightly more warped. The null hypotheses that the screened 
and unscreened distributions are identical is accepted only 
with a probability ^ 0.88. In Appendix Gl we show that this 
difference can be attributed to size or distance bias. In Ap- 
pendix G4 we show that warp distributions remain similar 
even if we expand the sample to include all types of observed 
warps. 

The distance-dependent effect, described in Appendix 
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Figure 7. An example of the estimation of warping from a disk galaxy image. The left panel shows the estimated warp curve. The red 
vertical line shows the x-center of the galaxy. The x-axis is given in the unit of half light radius (hr) of the galaxy. In the right panel we 
superimpose the warp curve on the optical image of the galaxy. 
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Figure 8. The distribution of warps for screened and unscreened galaxies for fjiQ = 2 x 10~^. The left panel shows the dimensionless 
warp estimated using Eq.ll and the right panel shows the estimate in kpc using Eq.l2. No significant difference is found between the 
distributions of warps of screened and unscreened galaxies in the left panel. A marginal difference is found on the right with more warp 
for unscreened galaxies. This discrepancy does not remain if we account for size bias as discussed in the text. 



CI, suggests that the ideal way to compare the W2 of 
screened and unscreened samples is to consider only galax- 
ies of similar physical size. The left panel of Figure 9 shows 
the warp W2 versus the half light radius of galaxies. As ex- 
pected, the warp is larger for bigger galaxies. Small points 
indicate measured values, and the large points show the av- 
erage and dispersion within 0.5kpc bins in half-light radius 
(see Appendix A). The lines shows the best fit linear rela- 
tion between the half light radius and the warp, with slopes 
listedin the top left corner. We see that when comparing 
galaxies of similar size, the warp is equivalent: this indicates 
that the difference of distributions in the right panel of 8 as 
due to this size bias. 



5.3 Comparison with theory 

Our analysis of galaxy warp shows that there exists no signif- 
icant difference between the occurrence of warps in screened 
and unscreened samples. Here we use this non-detection to 
place limits on strength of the fifth force. Given the low lev- 
els of observed warps, we assume that any contribution due 
to modified gravity will be added linearly to astrophysical 
effects which are estimated from the screened sample: under 
this assumption, to determine the warp associated with MG, 
we subtract the screened value from the unscreened value. 

We find the average warp {{wi)) and error on the aver- 
age using the method described in Appendix A. We find 
that (wi) = 0.035 ± 0.014 for unscreened galaxies and 
(wi) = 0.032 ±0.020 for screened galaxies (i.e. the warp due 
to astrophysical process is a 3% effect) for /rq = 2 x 10~^. 



© 0000 RAS, MNRAS 000, OOO-OOO 



12 Vikram, V. et al. 




Half light radius (kpc) Half light radius (l<pc) 



Figure 9. The left panel shows the strength of the warp in kpc as a function of galaxy size. The small black and red dots represent 
individual measurements of screened and unscreened galaxies for fjiQ = 2 x 10~^. The large circles with error bars show the average 
warp. The lines are the best fit linear relation to these averages. The slopes of the lines are shown in the top left corner. The right panel 
shows the estimated contribution to the warp due to MG, i.e., the difference between warp of unscreened (wu) and screened samples 
{ws). The blue and green lines shows the predictions with /S.G/G = 1 and AG/G = 1/3. Given the error bars, we are unable to test 
these models. 



This implies that the estimated warp in a modified gravity 
scenario (wmg) is as low as 0.003 ±0.024 (0.3%). See Table 
CI for the results of other background field values. We can 
use this in conjunction with a theoretical expectation of the 
warp strength to place loose constraints on modified grav- 
ity. This expected warp can be estimated as described by 
Jain & Vanderplas (2011); the result is sensitive to the mass 
profile of the galaxy, as well as the strength of the external 
potential gradient. Unfortunately, both these are difficult to 
determine from the data. 

To estimate the external potential gradient, we follow 
Jain & Vanderplas (2011) and assume that a galaxy will ac- 
quire a velocity of ^ 100 km/s in 3 Gyr. This gives an aver- 
age acceleration of ^ 10~^^ km/s^ We also need to consider 
the direction of this force at the galaxy to get the effective 
force to generate the warp. We marginalize over the angle by 
averaging 1000 random angles between -90 deg and 90 deg. 
The mass at a given radius of the galaxy is found by the 
assumption of cored isothermal sphere. Based on these as- 
sumptions we generate warp curves for these galaxies. This 
analysis leads to a theoretical prediction of the modified 
gravity warp contribution of wmg = 0.0031q qq^ (wmg — 
OmQtoml) for AG/G = 1/3 {AG/G = 1), with errors in- 
dicating the scatter in predicted values. 

The right panel of Figure 9 shows the average values 
of W2, as well as the theoretical predictions outlined above, 
as a function of galaxy size (these values are tabulated in 
Table CI). The blue and green lines show the expected warp 
for AG/G — 1 and 1/3, respectively. The large error in 
the measured values is due primarily to the fact that we 
are looking for a relatively small effect in a small sample of 
galaxies. The figure indicates that the measured warp is not 
precise enough to distinguish between modified gravity and 
GR, which corresponds to zero warp in this figure. 



5.4 Future Directions 

Our study of galaxy warps cannot rule-out either the predic- 
tions of MG or of GR. The two limitations are the relative 
scarcity of objects with kinematic data (like ALFALFA), 
and the uncertainty in halo parameters that goes into the 
warp prediction. 

Based on our estimates of uncertainties in the measured 
and predicted warps we can estimate the number of galax- 
ies with kinematic data required to reduce the error in the 
observed warp by a factor of 8 (so that the measurement 
error is comparable to the uncertainty in the predictions). 
We need ^ 8, 000 dwarf galaxies to test fno = 2 x 10~^ 
and ^ 20, 000 galaxies for /rq = 1 x 10~^ (calculated from 
third and fourth lines of Table CI). This assumes that the 
halo parameters have the same accuracy as we have in this 
paper. Note also that even though the number of galaxies 
required to test /rq = 2xl0~^isa factor of two less than 
that of /ro = 1 X 10~^, those galaxies should be very low 
mass and isolated, and therefore will require deep surveys. 
Another avenue of improvement is to obtain better optical 
images from current surveys (e.g Dark Energy Survey). This 
imaging should be supplemented by spectroscopic or 21cm 
surveys which can measure rotation velocities. 



6 STELLAR AND GASEOUS VELOCITY 
ROTATION CURVES 

In this section we explore the use of kinematic measures 
to detect effects of modified gravity. There are two classes 
of expected effects, both predicted by Jain & Vanderplas 
(2011). First, in the presence of a background potential gra- 
dient, a galaxy with a rotation axis perpendicular to the 
gradient will exhibit a distorted stellar rotation curve. Sec- 
ond, within the galaxy, the unscreened material will feel a 
stronger force and thus have its rotation curve enhanced by a 
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Figure 10. An example of rotation curves for two galaxies, mea- 
sured in Ha (dots with error bars) and HI (solid lines). The data 
are taken from Swaters et al. (2009). 



factor of y^l + AG/G. Both these measures depend on find- 
ing an accurate tracer of the kinematics of main sequence 
stars. We discussed Ha as a potential tracer in Section 3.2; 
we return to this subject here before proceeding with one of 
these tests. 



6.1 HI-Ha velocity difference 

Historically, Ha emission has been assumed to trace the stel- 
lar rotation curve. Because it originates in HH regions near 
hot stars, this assumption is well founded. However, as we 
argued in Section 3.2, this assumption may not hold in a 
modified gravity scenario. Because of the potentially com- 
plicated dynamical scenario of stars which lag behind the 
HH regions they create, however, simple arguments such as 
these may be incomplete. 

In this section, we address this question empirically by 
comparing the Ha rotation curve to that of HI, which is asso- 
ciated with unscreened gaseous component. Any similarity 
between these two within unscreened galaxies would indi- 
cate the safety of our assumption that Ha emission traces 
the unscreened galaxy component. 

In order to estimate the difference between HI and Ha 
rotation curves, we extract HI and Ha data from Figure 3 of 
Swaters et al. (2009), which contains 22 galaxies. The errors 
in vhi were obtained by private communication with Swa- 
ters. Each HI observation has one or two Ha counterparts 
(see the references therein); in addition, we have Ha data 
for 14 galaxies from GHASP. Note that all galaxies in the 
sample can safely be assumed to be unscreened due to their 
low mass and lack of massive neighbors. In Figure 10 we 
show the comparison of HI and Ha rotation curves for two 
galaxies taken from Swaters et al. (2009): the galaxy in the 
left panel has HI > Ha, while the galaxy in the right panel 
has HI < Ha. 

To quantify the difference between HI and Ha rotation 
curves statistically we first average both sides of the HI and 
Ha rotation curves separately about a common center to re- 
duce the observational uncertainty in the measured rotation 
curves, and then compare the HI and Ha curves to estimate 
the weighted average velocity difference, A^; = vhi — VHa- 
If both Ha and HI rotate similarly, then the /\v will be ap- 
proximately normally distributed about zero. In Figure 11 
we plot the distribution of vhi — VHa in black for Swaters 
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Figure 11. The difference between HI and Ho; rotation veloci- 
ties. We plot the distribution of A-u = vhi — VHa for the Swaters 
et al. (2009) data (black histogram). The red histogram includes 
individual errors in the Av measurements. We over plot in blue 
the best fit Gaussian. The mean value of the difference is indi- 
cated as well, showing that the rotation curves agree to better 
than 2 km/s. 



et al. (2009) data. When more than one Ha measurement 
is provided for the same galaxy, we weight equally the con- 
tribution of each Av to the mean. Accounting for errors in 
the observed values gives the red histogram.^ The best-fit 
Gaussian is indicated by the blue curve. 

We obtain (A^;) = 0.11 ± 1.50 db 5 km/s (statistical and 
systematic error) when we use HI and Ha from Figure 3 of 
Swaters et al. (2009). Combining different Ha observations 
and testing different ranges in the velocity curve leads to 
a 2 km/s contribution to the systematic error; the remain- 
ing 3 km/s comes from the potential for asymmetric drift 
& non-circular motion to affect the relationship (Swaters 
et al. 2009; de Blok & Bosma 2002). The errors additionally 
include those due to both measurement and non-circular 
motions, which we found by taking the difference between 
approaching and receding sides. See Swaters et al. (2009) for 
details. 

With these affects accounted for, we find that the HI 
and Ha rotation curves are consistent. This consistency con- 
firms that it is safe to assume that Ha traces the unscreened, 
gaseous component of the galaxy even in the presence of a 
fifth force. 



6.2 Asymmetry in Ha rotation curves 

In the previous section, we show that if MG is true, then 
Ha empirically traces the unscreened gaseous component of 
the galaxy. This implies that the Ha rotation curve will be 
nearly symmetric about the center of dark matter halo. On 
the other hand, the dark matter halo will be offset from 



^ In order to take into account individual errors and intrinsic 
dispersion to calculate the mean and the error we follow the pro- 
cedure described in Appendix A. 
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Figure 12. The Ha rotation curve asymmetry a = \Av\/v is 
shown as a function of rotation velocity for galaxies from the 
GHASP survey. Black (red) points indicate screened (unscreened) 
galaxies for fjiQ = 10~^. The mean and standard deviations of 
each distribution are shown by the large dots at the left of the 
plot. The shapes of individual points indicate the morphology 
of the rotation curves: diamonds indicate smoothly asymmetric 
curves, squares indicate barred galaxies, and empty circles indi- 
cate irregular galaxies. See text for details. 



the stellar component and therefore the Ha curve will be 
asymmetric with respect to the optical centroid, with an 
asymmetry on the order of 10-15km/s. We note that this 
asymmetry is a manifestation of the offset between the stel- 
lar and dark matter halo, and is distinct from the stellar 
rotation curve asymmetry predicted by Jain & Vanderplas 
(2011). 

In this section we use Ha rotation curves to test whether 
gaseous rotation curves are asymmetric about the stellar 
centroid, by comparing the approaching and receding sides 
of the Ha rotation curves about the associated stellar cen- 
troid. Any positive signal may indicate the presence of fifth 
force. 



6.2.1 Data and visual classification of rotation curves 

We use Ha rotation curves for nearly 200 disk galaxies from 
the Gassendi Ha survey of spirals (GHASP) archive (see 
Epinat et al. 2008, and references therein) to search for 
asymmetry in rotation curves. GHASP is an Ha kinematic 
survey of spiral and irregular galaxies. This survey has ob- 
served ^ 220 galaxies, and provides the rotation curves for 
both approaching and receding sides separately. We correct 
the velocities for inclination by assuming thin circular disk. 

We visually inspect these rotation curves and classify 
them into symmetric, smoothly asymmetric and irregular. 
We show a sample of symmetric curves in Figure D3 and 
smoothly asymmetric in Figure D4 in Appendix D. In a 
modified gravity scenario, we expect the gaseous rotation 
curves to be smoothly asymmetric about the stellar cen- 
troid (Jain & Vanderplas 2011). We remove all the irreg- 
ular rotation curves which have too few observed points 
and/or small-scale discontinuities: these are primarily as- 
sociated with small, dim and/or morphologically irregular 



galaxies, and the rotation curve shapes are more likely due 
to small-scale astrophysical effects than to a fifth-force type 
effect. We find that in both the screened and unscreened 
sample, symmetric curves outnumber smoothly asymmetric 
curves by a factor of four. We emphasize, however, that the 
galaxies with asymmetric Ha curves may be ideal for a fu- 
ture systematic search for the dynamical effects of modified 
gravity. 



6.2.2 Analysis 

We quantify the asymmetry in rotation curves for each 
galaxy using the abolute value of the weighted average of 
difference between approaching and receding sides of the 
rotation curves (|Av|) about the stellar centroids (Epinat 
et al. 2008). This asymmetry is better seen near the cen- 
ter of the galaxy and therefore we focus on the inner part of 
the rotation curves. However, we find that our results do not 
depend significantly on the location of the measurement. Fi- 
nally, we define the asymmetry coefficient a = \Av\/\vmax\ 
where Vmax is the maximum rotation velocity, since we ex- 
pect any asymmetries to scale with galaxy mass. The error 
in the asymmetry includes both measurement error and in- 
trinsic dispersion da, which can be same order of magnitude; 
we follow the method described in Appendix A. In a similar 
manner as the previous sections, we assume that the asym- 
metry of unscreened galaxy, is a linear combination of 
, the modified gravity contribution, and a^, the contri- 
bution from standard astrophysical forces. 



6.2.3 Results 

We summarize the results in Figure 12. We plot a as a 
function of maximum rotation velocity of the galaxies for 
fno = 10~^. There are 68 unscreened galaxies (shown in red) 
and 58 screened galaxies (shown in black). We compute the 
asymmetry within the central region of the rotation curve, 
up to 30% the optical radius of the galaxy. 

The shapes of the individual symbols indicate the 
visually-determined galaxy morphology: the diamond sym- 
bols indicate galaxies with visual asymmetry in the rotation 
curve which look interesting under modified gravity theo- 
ries. Square symbols indicate barred galaxies and empty cir- 
cles represent irregular cases. There is no clear correlation 
between these morphologies and the measured asymmetry. 
The filled points on the left of the plot show the mean and 
standard deviation of the asymmetries for the screened and 
unscreened components. We find that the asymmetry for 
both screened and unscreened galaxies is around 7% which 
corresponds to asymmetries in the range of 10-20 km/s. Sub- 
tracting these mean values implies an asymmetry produced 
by modified gravity of a^^ = -0.001 ± 0.012. If we in- 
stead assume fno = 10 (24 unscreened galaxies and 102 
screened galaxies) we find a^^ = 0.024 ± 0.015. In either 
case, a 10% asymmetry predicted for modified gravity can 
be ruled out at approximately the 5-cr level. 

Note that these results do not attempt to correct for 
asymmetric drift or non-circular motion. We have further 
analyzed this data and found no correlation between asym- 
metry and the number of neighbors or distance to the near- 
est neighbor, which would be present if the asymmetries 
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were caused by a fifth force. Garrido et al. (2005) also stud- 
ied asymmetry in GHASP Ha rotational curves. They con- 
clude that both isolated and softly interacting galaxies have 
similar level of asymmetry, which might be due to inhomo- 
geneous massive star formation regions. 



7 DISCUSSION AND CONCLUSIONS 

We have explored several astrophysical tests of modified 
gravity which rely on partially-screened dwarf galaxies in 
which fifth force effects act on diffuse components (gas, dark 
matter, and giant stars), while more compact components 
(main sequence stars) are screened from these effects. These 
tests probe regions of the parameter space of chameleon type 
gravity theories which have have not been probed by cos- 
mological tests or laboratory experiments. The observables 
include offsets between the stellar and gaseous components, 
warping of stellar disks, and asymmetry in the gaseous ro- 
tation curves. Our study relies on a variety of previously- 
published observational data. 

7.1 Summary of results 

• Offset between HI gas and stars: The offset study based 
on optical and HI gas centroids shows that it is not possible 
to distinguish between gravity theories using the current HI 
data. The accuracy of the astrometry is principally limited 
by the large beam size of the 21cm observations. However, 
with a larger catalog of candidate galaxies and/or more ac- 
curate radio astrometry, this test may provide competitive 
constraints in the future. In this section we summarize our 
results, discuss systematic errors and examine prospects for 
improved tests. 

• Offset between optical and dynamical center: An alter- 
native means of measuring this expected offset, the differ- 
ence between the optical centroid and the kinematical cen- 
troid of Ha, is much more precise and limits the component 
of modified gravity scatter to ^ 2% of the optical radius of 
the galaxy, ruling out 5% scatter at the 3-cr level. 

• Offset between red giants and main sequence stars: We 
show that the offset between the partially-screened red giant 
population and the fully- screened main sequence population 
is too sensitive by systematic effects to provide any meaning- 
ful power to distinguish between gravity theories: the data 
is consistent with both GR and MG at a 2-cr level. This test 
has the potential to improve as more stellar populations are 
observed within nearby dwarf galaxies. 

• Warps in edge- on galaxies: The warp study relies on 
measurements of the optical morphology of stellar disks of 
late-type galaxies observed by SDSS. We find that the warp 
measures for this galaxy sample are consistent with both 
GR and modified gravity theories: the main limitation is 
the relatively small number of suitable galaxies for which 
kinematic data is available in order to estimate the screening 
level of the galaxies in the sample. 

• Rotation curve tests: In the final test we used the 
GHASP galaxy sample to estimate the asymmetry in the 
Ha rotation curves, which arises due to the offset between 
the screened stellar component and the unscreened gaseous 
component. We find that the asymmetry attributable to 
modified gravity is as low as (-0.1 =b 1)%, compared to the 



expected value of 10% under favorable conditions. This im- 
plies that the modified gravity models which predicts 10% 
asymmetry can be rule out at 5cr level. However, this test 
carries two major uncertainties: systematic uncertainties in 
the data due to effects of asymmetric drift and non-circular 
motions, and systematic uncertainties in the theoretical ex- 
pectation due to insufficient knowledge of the internal halo 
parameters and external environment of each galaxy. 

Finally, we examined individual rotation curves from 
the GHASP galaxy sample and classified them based on 
the different signatures in the curves (see Appendix D). 
We found that the majority of these galaxies have fairly 
symmetric curves regardless of their local density /potential. 
However, we identified a few galaxies with interesting fea- 
tures in the rotation curves. Some of these galaxies have 
other potential fifth force effects like U-type warps and op- 
tical offset from that of HI. More detailed studies of these 
individual cases can shed light on whether the features may 
be attributable to fifth-force effects or more traditional as- 
trophysical processes. 

7.2 Systematic errors and caveats 

The tests performed here involve aggregate statistical com- 
parisons between the unscreened galaxy sample and the 
screened (control) sample. One of the major assumptions 
in this type of study is that, apart from screening level, the 
two samples are statistically similar. As noted above, for 
many of the tests this division correlates with galaxy size. 
In some cases we attempt to minimize any potential bias by 
looking only at galaxies near the self-screening threshold (at 
the cost of weakening the signal). However, a more precise 
result could be obtained given a large sample of low-mass, 
environmentally screened galaxies. Galaxies in the two sam- 
ples would then have the same distribution in size and mass. 
So we recommend that future studies in this direction prior- 
itize new observations of low-mass galaxies in both the field 
and near large neighbors. 

Another important caveat is that we use a large, het- 
erogenous data set to perform the above tests. Because the 
data were gathered without this type of test in mind, there 
was little specific attempt to control systematic errors that 
can affect explorations of modified gravity. One extreme ex- 
ample of this is the constraint based on separation of HI and 
optical centroids: the strength of the test is sensitive to the 
characteristics of the error in centroid determination: even a 
10% overestimate of the error can impact the resulting con- 
straints. Observing campaigns specifically designed for this 
sort of modified gravity test could yield homogeneous data 
sets with such systematic effects minimized. 

The third major caveat comes from uncertainties in our 
estimate of the external potential at the locations of dwarf 
galaxies. We have estimated the Newtonian potential as well 
as the fifth force, estimated using the masses and locations 
of neighboring galaxies. This estimate is subject to the un- 
certainties in mass estimation of galaxies. It can be improved 
in a number of ways: completeness over the neighborhood of 
the dwarf galaxy sample, mass estimates using spectroscopic 
or 21cm observations, and better modeling of the mass and 
scalar field distribution. The improved modeling can be done 
by explicitly solving for the scalar field at each location, in- 
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eluding the presenee of underdense regions that also con- 
tribute to the fifth force. For some of the tests only a ~ 10 
Mpc sized region needs to studied in detail, so the observa- 
tions and modeling are tractable. 

7.3 Future prospects 

Here we discuss specific ways to sharpen each of the tests 
we have carried out and summarized above in subsection 
7.1. Before we get to these tests, we note that a variety of 
experiments can be used to test gravity theories, ranging 
from laboratory to cosmological scales. For chameleon type 
scalar-tensor theories, the first in our series of papers (Jain 
et al. 2012) used cepheid variable stars and other distance 
indicators to place stringent limits on the two parameters 
of the theories: the coupling a (related to AG/G in this 
paper) and the background field value (equivalent to fno 
in this paper). The upper limit for f{R) models obtained 
using cepheids correspond to AG/G — 1/3; fno ^ 5 x 10~^, 
or a fifth force range below 1 Mpc at cosmic mean density: 
these are far more stringent than limits from large-scale, 
cosmological tests as discussed by Jain et al. (2012). For 
tests using other tracers in the Milky Way, see Brax & Davis 
(2013) and reference therein. 

Dwarf galaxy tests currently probe larger coupling val- 
ues. However the field values probed can be lower, especially 
with the fainter dwarf galaxies within 10 Mpc of the Milky 
Way. The limits on both parameters can be improved signif- 
icantly for all the tests we have considered. A brief account 
for each test follows, though detailed studies are needed to 
plan future observations. 

• Offset between HI gas and stars: For this test to pro- 
duce competitive constraints, it is crucial to improve the 
precision of the centroid measurement, either through tele- 
scopes with a smaller beam size, or through the stacking of 
many more galaxy images. Using the ALFALFA data, the 
error in the offset estimate is ~ 1 kpc; to get an error of 100 
pc we need to stack ~ 10^ galaxies. Furthermore, the char- 
acteristics of the error in the centroid must be understood 
at the few percent level. A better way forward would be to 
get improved resolution radio maps (e.g. via VLA observa- 
tions) , which would then relax the required sample size and 
lower the systematic uncertainty. 

• Offset between optical and dynamical center: increasing 
the number of Ho; rotation curves to a few hundred low- 
mass galaxies would be sufficient to reduce the uncertainty of 
this measurement. More important is sophisticated model- 
ing of the environment of each galaxy to determine screening 
levels and the strength of the expected signal. 

• Offset between red giants and main sequence stars: Cur- 
rently, the red giant displacement test is dominated primar- 
ily by systematic uncertainties. As the signal in this test is 
weaker than the other probes, the best approach would be 
to focus on a limited number of galaxies with a very well 
known external force (i.e. well-observed and well-modeled 
local environment). The uncertainty in the external force is 
the limiting factor if one wants to constrain modified grav- 
ity. Another potential approach would be to use integral- 
field galaxy spectra and estimate the contributions of red 
giant and main sequence populations as a function of loca- 
tion in the galaxy. A difference in these centroids could be 



attributed to modified gravity, using control tests similar to 
those in this paper. This test is applied to a small number 
of very small but nearby galaxies (closer than 10 Mpc), so 
detailed studies of individual galaxies are feasible. 

• Warps in edge- on galaxies: To reduce the observed er- 
ror in the warp measurements to the level of uncertainty in 
the predicted values requires high-quality images of about 
10,000 dwarf galaxies in very low density environments 
and 20,000 dwarf galaxies in low to moderately dense en- 
vironments for f{R) theories with fno = 2 x 10~^ and 
fno = 1 X 10~^, respectively. These numbers are feasible 
with current and ongoing galaxy surveys. With this increase 
in data, however, the limiting factor will be the knowledge 
of the distribution of halo parameters within galaxies of this 
size: the predicted strength of the observed warp is highly 
dependent on this mass profile. Rotation curves that extend 
to the inner parts of the disk can constrain the mass profile. 

• Rotation curve tests: We have examined asymmetries 
in rotation curve data. In addition, a direct measurement 
of AG/G can be made by comparing HI rotation curves 
with stellar rotation curves (section 3.1). However the stel- 
lar rotation curve must be measured using stellar absorp- 
tion lines which arise very close to the surface of dwarf stars 
(not via Ha as is common practice; see discussion in sec- 
tion 3.2). These curves are difficult to observe: depending 
on brightness, each galaxy may require several hours of ob- 
servations with large telescopes. This test has the advantage 
that it does not depend on a detailed dynamical model or 
on the precise characteristics of the environment of the un- 
screened galaxy: for this reason, high-precision observations 
of a handful of the unscreened galaxies can lead to strong 
constraints on modified gravity theories, free of many of the 
systematic uncertainties discussed above. 
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Figure Bl. Color-magnitude diagram for galaxy DD0187, using 
filters F606W and F814W. The red region shows the selected red 
giants, that reside at the bright tip of the red giant branch. 



APPENDIX A: ALGORITHM TO FIND 
INTRINSIC DISPERSION 

In most of the probes that we consider in this paper, there 
is an intrinsic dispersion associated with the estimated pa- 
rameter due to astrophysical systematic effects that can be 
larger than the error in individual measurements. Because 
we weight the data by the inverse variance while calculating 
the mean value, including an intrinsic dispersion keeps data 
points with small measurement errors from being inappro- 
priately over- weighted. It leads to an improved estimate of 
the mean and the error in the mean value. Here we describe 
our approximate procedure to estimate the intrinsic disper- 
sion, and the improved mean and the error in the mean. 

We have set of N observed values xi and errors cr^. We 
calculate the weighted average as 



Y^i XjWj 



(Al) 



where Wi — 1/crf. Now, we add a constant systematic error 



(7s to all the individual errors as (7^^^^ 



ai + as which 



accounts for the intrinsic dispersion. For a Gaussian random 
variable this satisfies: 

{xi - xY 



N 



(A2) 



Now we recalculate the mean x using Eq. Al with Wi — 
l/^efF,i- This process continues until it converges. Finally, 
the error on the average x will be calculated by 

1 



2 



(A3) 



APPENDIX B: RED GIANT STARS 

In this section we summarise the potential systematic errors 
related to test based on red giant stars. In Figure Bl we 
show an example of color-magnitude diagram to illustrate 
the selection of red giant stars. 



In the left panel of Figure B2 we explore whether the 
observed displacement between the red giant and main se- 
quence populations is correlated with the distance to the 
galaxy. The data shows a small positive correlation between 
distance and displacement. To check whether this might af- 
fect the conclusions of Section 4.3, we divide the sample in 
different distance bins and reanalyze them in individual bins. 
We find that our conclusions do not change significantly even 
if we consider galaxies at similar distances. 

In the right panel of Figure B2 we plot the correlation 
between the displacement and a measurement of asymme- 
try. The asymmetry is defined as the difference in centroids 
of main sequence stars at two different isocontours of stel- 
lar number density describing inner and outer part of the 
galaxy. We see that red giant displacement is strongly cor- 
related with this asymmetry measurement, indicating that 
the effect may be related to the internal structure of the 
galaxy rather than being due entirely to a fifth-force effect. 
We note, however, that similar morphological asymmetry 
may occur in the presence of modified gravity as well. How- 
ever, if we find the presence of similar effects in a screened 
galaxy sample then we will able to separate the modified 
gravity effects from internal structure of the galaxy. 

We also checked if the displacement is correlated with 
the number of stars (i.e. shot-noise), the ratio of red giants 
compared to the number of main sequence stars, the lumi- 
nosity of the galaxy (indicative of mass and size), the HI 
line width (also indicator of mass), filters used to determine 
red giants or magnitude limit in the whole sample which de- 
pends on the camera used. None of these observables show 
appreciable correlation with the displacement of red giants. 

We examined the effect of differences in the area that we 
used to calculate the centroid of main sequence stars and red 
giants. Usually the galaxy is not well centered in the image, 
that is the reason why we have to study systematics in the 
definition of the galaxy limiting area. Although it is easier 
to define isocontours in number counts using all the stars, 
there are a lot of interlopers (faint stars) that do not belong 
to the galaxy. We recalculated the centroid after gradually 
removing the faintest stars. We include galaxies only if their 
centroids converge after removing 20% of faintest stars. We 
also test if crowded regions can bias the result based on 
flags defined in Dalcanton et al. (2009). Due to longer ex- 
posure times, large number of faint foreground star will be 
detected which leads to poor photometry of individual stars. 
Moreover, we have calculated the centroid of main sequence 
stars weighting by luminosity, and still the displacement of 
red giants is consistent. Finally, we change our definition of 
red giants by going 0.5 magnitudes fainter, and we exclude 
bright stars from the definition of main sequence stars, and 
conclusions remain the same. 

We also checked if the displacement is related to the 
number of stars (shot-noise), the ratio of red giants com- 
pared to the number of main sequence stars, the luminosity 
of the galaxy (indicative of mass and size), the HI line width 
(also indicator of mass), filters used to determine red giants 
or magnitude limit in the whole sample which depends on 
the camera used. None of these observables are correlated 
with the displacement of red giants. 
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Figure B2. We plot distance vs displacement in the left panel to study systematic effects due to lower resolution for more distant 
galaxies. In the right panel, we plot the correlation between red giants displacement and degree of asymmetry of the galaxy (the straight 
line shows y=x). Note how the asymmetry is bigger than the displacement in almost all the cases, suggesting that the displacement is 
just a result of internal structure. 



APPENDIX C: MEASUREMENT OF WARPS 

In this section we describe the systematics effects related to 
the measurement of warps in section 5. We explore the bias 
due to distance effects, inclination effects, and environment. 
Finally, we propose a third warp estimator and show that 
the conclusions of Section 5 are unchanged. 



CI Distance Bias 

A number of systematic effects may be involved in the warp 
analysis. The first is bias related to distance. In our sam- 
ple there exists a strong correlation between distance and 
average physical galaxy size: this is due to the fact that it 
is difficult to observe physically smaller galaxies at higher 
distances. If the measured warp depends on the size of 
the galaxy, then the comparison between screened and un- 
screened samples should be done as a function of distance 
to limit systematic effects due to the different samples. We 
check the distance bias by plotting the warp parameter 
against distance to the galaxies. The left panel of Figure 
CI shows that the normalized warp parameter, wi, is inde- 
pendent of the distance to the galaxy, as expected. 

On the other hand, the estimator W2 does depend on 
the size of the galaxy: two galaxies with identical wi curves 
will differ in W2 , depending on the size of the galaxy. There- 
fore, the size-distance relation implies that the estimated 
values of W2 will be higher at larger distances. The right 
panel of Figure CI shows that this is indeed true. It should 
be noted that we have a systematically higher fraction of 
screened galaxies at smaller distance. We found that these 
galaxies belong to the Virgo cluster at distance ~ 17 Mpc. 
Therefore, when we consider all the galaxies irrespective of 
their size, the screened sample has a smaller average value of 
W2 compared to that of unscreened galaxies due to the large 
fraction of screened galaxies at smaller distances. That is, 
the marginal difference between W2 distributions of screened 
and unscreened samples can be explained by this effect (com- 
pare figure 8 and Figure CI). 




0.3 0.4 
Axis ratio 

Figure C2. The distribution of axis ratio for screened and un- 
screened galaxies with fjio = 2 x 10~^. We find no difference 
between the two samples. 



C2 Inclination Bias 

Another potential systematic effect may be related to the in- 
clination of the galaxy. While the warp can be well-estimated 
in edge-on galaxies, the apparent warp may be nearly zero 
for warped face-on galaxies. This is the source of the bias 
introduced by inclination. If the unscreened galaxies are sys- 
tematically more inclined than the screened galaxies, then 
the warp parameter estimation will be more robust for un- 
screened galaxies. We therefore, check whether there is any 
systematic difference between the distributions of axis ratio 
for screened and unscreened galaxies. Figure C2 shows that 
the axis ratio of both screened and unscreened galaxies are 
distributed similarly. So we conclude that inclination affects 
the screened and unscreened galaxies equally, and therefore 
is unlikely to bias our result. 
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Figure CI. Systematic effects involved in the warp test. The left panel shows the warp (wi) vs distance (Eq. 11), right panel shows 
warp in kpc vs distance (Eq. 12) for screened and unscreened galaxies. If the screened sample is systematically closer than unscreened 
galaxies, then we measure a lower value of warp for screened sample compared to unscreened sample, which goes in the same direction 
as the modified gravity prediction. 



C3 Environmental Bias 

An additional test we can perform with the data is to check 
whether there is any noticeable difference in the fraction 
of galaxies displaying U-shaped warps between the screened 
and unscreened samples. In the modified gravity scenario we 
expect that the fraction of U-shaped warps will be larger in 
the unscreened sample than in the screened sample. For this 
test we classify galaxies as U-shaped only if the galaxy has 
more than 3% warp, i.e. if wi > 0.03. As above, if the signs 
of left and right side warp are different then we classify them 
as U-type. 

Figure C3 shows that the fraction of U-shaped warps 
does not depend on the screening level. This implies that 
any physical process which causes the formation of U-shaped 
warp works irrespective of the local density of the galaxy. 
It should be noted that this result does not vary with the 
choice of lower threshold value of the warp (i.e. wi > 3%) 
to define the U-warp. In the following section, we further 
show that the fraction of S-shaped warps is also consistent 
between the screened and unscreened sample. 
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Figure C3. The fraction of galaxies displaying U-shaped warps 
vs the magnitude of the external potential. In this plot, galaxies 
are classified as U-shaped only if the warp is greater than or equal 
to 3%. 



C4 Distribution of warps 

In Section 5 we have shown that the distribution of U-type 
warp does not depend on the the environment. However, we 
include only U-type warp in that analysis. Here we check 
whether the conclusion in section 5 still holds if we include 
all kind of warp in the analysis. For this purpose we de- 
fine the effective total warp (we) of galaxies as the abso- 
lute value of the difference between the warp of left and 
right hand sides of the galaxy. As we explained before the 
U-type galaxies will have opposite signs for left and right 
side waprnesses and S-type/N-type will have same signs for 
both sides. Therefore, our definition of We will leads to larger 
value for a symmetric U-type warp and zero for other cases. 
However, the other type of warps can also give smaller non- 
zero values for We depending on the asymmetry of the warp 
curve. The same applies to U-type warp where the We can 
be smaller for weak warp curves. We find that almost all 
galaxies, except U-type warps, have We < 0.1 with warp 



definition Eq 11. It is shown in the left panel of Figure C4. 
Also, we found that the distribution of We does not have 
any dependence on the environment and it is shown in the 
middle panel of Figure C4. We also found that there is no no- 
ticeable difference between the distributions of S-type warps 
in screened and unscreened regions and it is shown in the 
right panel of the figure. We found similar results if we use 
Eq 12 to define warp. In that case We = 0.5 kpc separate 
U-type warps from the rest. 

In Figure C3 we show the fraction of galaxies with U- 
type warps in different environments. We check whether the 
fraction of normal galaxies and fraction of S-type galaxies 
have any dependence on their environment and the result is 
shown in Figure C5. We classify a galaxy as normal if the 
warp is less than 5%. We find that the fractions of different 
types of galaxies is independent of their environment. The 
result is insensitive to the definition of normal galaxy. 
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Figure C4. Left : The cumulative histrogram of We (see text for the definition). The U-type warps can be distinguished with We > 0.1. 
Middle panel shows the distributions of We in screened and unscreened regions for fpiQ = 2x 10~^. The right panel shows the distributions 
of S-type warps in screened and unscreened regions with /^q = 2 x 10~^. 



Type fpiQ Number (^i) '^mg Prediction 

Ns Nu ws wu AG/G = 1/3,AG/G = 1 

All 2e-07 59 158 0.051 ± 0.023 0.050 ± 0.014 -0.001 ± 0.027 

All le-06 143 352 0.038 ± 0.014 0.046 ± 0.009 0.008 ± 0.017 

U 2e-07 48 78 0.031 ± 0.018 0.036 ± 0.015 0.006 ± 0.023 0.003g;0gf , 0.006g;gg^ 

U le-06 89 229 0.024 ± 0.012 0.029 ± 0.009 0.005 ± 0.015 0.0020;00f, 0.004g;gg| 
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Table CI. Summary of warp tests. The columns are: 1. Type: the type of warp included in the analysis. All - included all galaxies, U 
- used only U-type warps. 2. fjiQ - the value of fpiQ parameter used to classify screened and unscreened galaxies 3. Number: - A^^, Nu 
are the number of galaxies in screened and unscreened regions. 4. ws^ wu are the average value of the two wi or W2 in screened and 
unscreened galaxies 5. wmg - estimated warp due to modified gravity. 6. Prediction: the theoretical predictions for unscreened galaxies 
with AG/G = 1/3 and AG/G = 1. The upper part of the table is for the dimensionless warp parameter wi (defined in Eqn. 11) and 
lower part is for the parameter W2 (defined in Eqn. 12). In the second part of the table we group U-type galaxies based on their half 
light radius (hr) to compare the warp. See text for details. The estimated value of wmg is consistent with GR in all cases, but the error 
bars are too large to detect the MG predictions shown in the last column. 



APPENDIX D: EXAMPLES OF INTERESTING 
OBSERVATIONS 

We present some examples of optical images (Figures Dl & 
D2) and rotation curves (Figures D3 & D4) of unscreened 
galaxies from the literature which look interesting in the 
context of modified gravity. 
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Figure Dl. Left panel: Optical image of dwarf galaxy UGC 3740 which is unscreened for fjiQ = 10~^. The galaxy is asymmetric in the 
optical. The Ho; rotation curve also shows strong asymmetry. Right panel: Comparison of HI and optical images of UGC 4325. The HI 
image shows some asymmetry. However, this galaxy does not have any nearby neighbors to which we could attribute an external force. 
The images are taken from Begum et al. (2008); Swaters et al. (2002) 




Figure D2. Optical images of galaxies with the most asymmetric Ho; velocity curves: Left panel: UGC 7831 which is unscreened for 
fuQ = 10~^ and shows some warp. Middle panel: UGC 5721. Right panel: UGC 1913. The images are taken from Begum et al. (2008) 




Normal 



Figure C5. The fraction of normal, S-type and U-type galaxies 
in screened and unscreened regions. We find no noticeable dif- 
ference between the fraction of different types in screened and 
unscreened regions. The error in the plots are calculated by as- 
suming Poissonian noise on the number of galaxies. 
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Figure D3. Examples of Ho; symmetric rotation curves for GHASP galaxies. The x-axis is in kpc, and the y-axis shows the measure 
rotation velocity in km/s (before correcting for inclination). The left side points of the rotation curves are shown as red points for 
comparison to the right side. We indicate the name of the galaxy, Av defined in the text and the external potential due to neighbors 
within 3 Mpc. The potential is in red if the galaxy is unscreened, and black if it is environmentally screened, for fjiQ = 10~^. We label the 
inclination and the averaged maximum velocity of the galaxy once corrected for inclination. Diamonds show HI curves, when available. 
The blue vertical lines show the range used to calculate Av, based on its half-light radius. Note that for UGC 4325 the velocity in HI is 
larger than in Ha. 
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Figure D4. Same as Figure D3 for smoothly asymmetric curves, selected visually. These slow varying asymmetries could be created by 
modified gravity. 
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